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We estimate the effect of the superconducting gap anisotropy in the dispersive gap mode of 
phonons, which is observed by the neutron scattering on borocarbide superconductors. We nu- 
merically analyze the phonon spectrum considering the electron-phonon coupling, and examine 
contributions coming from the gap suppression and the sign change of the pairing function on 
the Fermi surface. When the sign of the pairing function is changed by the nesting translation, 
the gap mode does not appear. We also discuss the suppression of the phonon softening of the 
Kohn anomaly due to the onset of superconductivity. We demonstrate that observation of the 
gap dispersive mode is useful for sorting out the underlying superconducting pairing function. 

KEYWORDS: gap mode in phonons, anisotropic superconducting gap, electron-phonon coupling, borocarbide. 



1. Introduction 

In the materials with strong electron-phonon interac- 
tions, we can obtain the information of the electronic 
state by investigating phonon properties. The Kohn 
anomaly of the phonon dispersion comes from the nest- 
ing features of the band structure in simple metals via 
electron-phonon interaction. 1 ' It is known that the gap 
mode of phonons appears at w = 2A in the supercon- 
ducting state with a gap A. In 2H-NbSe2, the Raman 
scattering 2 ' observes the gap mode around zero phonon 
wave vector q = 0. The gap mode at q 7^ is reported 
by the neutron scattering on YNi 2 B 2 C. 3 ~ 5 ' There, the 
phonon dispersion shows the Kohn anomaly at the nest- 
ing wave vector Q ~ 7r(0.55, 0, 0), 6 ' and the new peak 
of the dispersive gap mode appears below the original 
phonon dispersion around q = Q. The neutron scatter- 
ing on LUM2B2C reports that the peak energy at the 
phonon dispersion is shifted to lower energy and that 
the peak intensity is enhanced below the superconduct- 
ing transition temperature T c . 7 ' Theoretically, the gap 
mode is obtained by considering the contribution of the 
electron-phonon coupling in the phonon sclfcncrgy. 8-11 ' 
However, previous theoretical works have investigated 
in the isotropic s-wave superconductivity case, and they 
have not considered the dispersive behavior of the gap 
mode. 

Band calculations suggest that borocarbide supercon- 
ductors are conventional BCS superconductors with the 
pairing interaction mediated by the electron-phonon in- 
teraction. 12 ~ 15 ) However, the direct observation of the 
superconducting gap of YM2B2C by photoemission spec- 
troscopy reports that the gap has large anisotropy and 
that the minimum of the gap seems to be zero as if it 
has point or line nodes. 16 ' And the impurity effects 
suggest that they are extreme anisotropic s-wave su- 
perconductors. This gap structure is also supported by 
the specific heat measurement. 17 ' 18 ' Thermal conductiv- 
ity experiment by rotating magnetic field suggests that 
the superconducting gap has point nodes. 19 ' Therefore, 
it is important to study how the superconducting gap 
anisotropy affects the dispersive gap mode of phonons. 
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In this paper, we investigate the effect of the gap 
anisotropy on the gap mode of phonons by calculating 
the dispersion and intensity of the gap mode around 
q = Q. As the gap anisotropy effect, we examine the 
case when the superconducting gap is suppressed around 
the node points, and the case when the order parame- 
ter changes the sign by the nesting translation on the 
Fermi surface as in the d-wave pairing. The phonon 
properties are studied by varying 2A from the weaker 
case to the stronger case relative to the phonon fre- 
quency And we discuss the experimental data observed 
in YNi 2 B 2 C and LuNi 2 B 2 C. To obtain the renormalized 
phonon Green's function in the random phase approxi- 
mation (RPA), we calculate the dynamical charge sus- 
ceptibility in anisotropic superconductors. In the cal- 
culation, we use the electron dispersion characterized by 
the nesting vector Q ~ 7r(0.55, 0, 0), which reproduce the 
Kohn anomaly. When the Kohn anomaly is strong, the 
phonon is softened at q ~ Q 1 indicating the transition 
to the charge density wave (CDW) state. This softening 
is suppressed by the onset of superconductivity. 

After explaining our formulation in §2, we present our 
results of the phonon spectral function and discuss the 
effect of anisotropic superconductivity in §3. The sup- 
pression of the phonon softening is discussed in §4. The 
last section is devoted to summary and discussions. 

2. Formulation 

As the Fermi surface shape of borocarbide supercon- 
ductors is complicated, 6 ' 15 ' we simplify the dominant 
large Fermi surface to the simple two-dimensional Fermi 
surface which reproduces the nesting property of nest- 
ing vector Q = (Q X7 0,0) ~ 7r(0.55, 0, 0) as in the boro- 
carbides. Then, we use the following dispersion of the 
electrons, 



K, 



-2i(cos k x + cos k y ) — At' cos k x cos k y 
-2t"(cos 2k x + cos 2k y ) — (i, 



(1) 



with the transfer between nearest (t), second nearest (t') 
and third nearest (t") neighbor sites in two-dimensional 
square lattice of the basal plane in tetragonal crystal. 
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The lattice constant is taken to unity. Here, we set 
t = 1, t' = —0.6, t" = 0.1 and the chemical potential 
^ = —1.05. The energy and temperature are scaled by 
t throughout this paper. The Fermi surface obtained by 
K k = is shown in Fig. 1 (a) . Since it has flat surface at 
k x ~ ±0.2757r, it gives good nesting for Q ~ 7r(0.55, 0, 0). 
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Fig. 1. (a) Fermi surface for t = 1, t' = -0.6, t" = 0.1 and 
fi = —1.05 in momentum region < k x < n and < k y < n. (b) 
The gap amplitude of the pairing functions along the Fermi sur- 
face. We plot \<f> s (k)\, |0d x2 _ B2 O)l = \4>s x2 _ y2 (fc)|, |</>ci X3 ,(fe)l = 

I'/'sxy COli an d I'/'pa; CO I 38 a function of 0, which is the angle cen- 
tered at the X point as shown in (a). The lower panels schemat- 
ically present the sign of the pairing function along the Fermi 
surface for the d x i _ y 2-wave (c), the d xy -wave (d) and the p x - 
wave (e) cases in the region — tt < k x < tt and — tt < k y < tt. 
The solid circles on the Fermi surface present the node points of 
the pairing functions. Arrows show the nesting vector Q. 

As for the pairing function (j)(k), we use <p s {k) = 1 
for isotropic s-wave pairing, 4>d x2 _ y2 (&) = cos k x — cos k y 
for d x 2_ !/ 2-wave, and (^^(k) = 2 sin k x sin k y for d xy - 
wave. The gap amplitude along the Fermi surface is 
shown in Fig. 1(b). The flat part of the Fermi sur- 
face related to the nesting is ranged for < 9 < 30°. At 
6 = 0, \4>d x2 _ y2 (k)\ has large gap, and \4>d xy (k)\ = due 
to the node structure. With increasing 9, \<pd us2 _ 2 {k)\ 
{\<pd xy (k)\) is decreased (increased). These gap depen- 
dences on 9 are averaged along the flat Fermi surface 
when they contribute to the phonon properties. 

It is also important to see the sign change of the pair- 
ing function through the nesting translation Q. In the 
d;j.2_ y 2-wave case, there are nodes along the diagonal 
lines k y = ik x . In the c^y-wave case, nodes are located 
along the horizontal lines k y — 0, ±n and the vertical 
lines k x = 0, ±n. The sign of 4>d x2 _ y2 {k) is the same at 
the both ends connected by the nesting Q, as schemati- 
cally shown in Fig. 1(c). But, the sign of <j>d (k) changes 
by the nesting Q, as shown in Fig. 1(d). 

For extreme anisotropic s-wave superconductors, we 
use <j> Sx2 _ 2 (k) = \0d !c2 _ y2 {k)\ and <j> Sxy (k) = \(pd xy {k)\. 
We call them as the s x 2_ y 2 case and the s xy case, re- 
spectively. They have the same superconducting gap as 
those of the d x 2 _ y 2 -wave and the e^-wave cases. But 
they do not change the sign of the pairing function un- 
der 7r/2-rotation. Then, the pairing function does not 



change the sign by the nesting Q. We also consider the 
p^-wave pairing case 4> P:i ,(k) = sinfc^, while it is not for 
the pairing of borocarbide superconductors. The pairing 
function Px (k) has nodes along the vertical lines k x — 0, 
±7r. As shown in Fig. 1(b), 4> Vx (k) has large gap at the 
flat Fermi surface < 9 < 30°. But, <p Plc {k) changes the 
sign by the nesting Q as shown in Fig. 1(c). Then, we 
can clarify the effect of the sign change by considering 
<t> Px {k). 

First, we calculate the dynamical charge susceptibility 
Xnn(q,u) in the superconducting state with the super- 
conducting gap Afc = A(T)(j)(k). This is given by 



1 
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K k+q K k - Rc{A fc+q A* } 

Ek+qEk 



x n(E k+q )-.f{E k ) | f(E k+q ) f(E k ) 
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f{E k+q ) 



E k + q E k 

f(E k )-l f(E k+q ) + f(E k )-l 



with E k = (K k 



E k+q + E k - uj E k+q + E k + u 

2 'Afcl 2 ) 1 / 2 and the Fermi distribution 
function f(E). We perform fc-summation within a Bril- 
louin zone. The temperature dependence of A(T) is de- 
termined from the BCS relation, 



1 - 2/Cgfc) 
2E k 



(3) 



where the pairing interaction V s (< 0) is determined by a 
given A(T = 0). 

The neutron-scattering profile of a phonon is approx- 
imately the phonon spectral function S(q 1 tu) 1 i.e., the 
imaginary part of the phonon Green's function, 



S(q,w) = -TmD{q,w) 



(4) 



The renormalized phonon Green's function D{q,w) is 
given by 

D^u)- 1 =D (q,oj)- 1 -U(q,u)). (5) 
In the phonon sclfcncrgy, polarizability is written as 

I%,w) = -\g q \ 2 Xnn(q,u) (6) 
in RPA. The bare phonon Green's function 

2uj n 



Do(q) = 



UJ 2 — LU q 



(7) 



The electron-phonon coupling constant g q and the 
phonon dispersion u q are treated as constants, since we 
consider narrow qr-region around Q. 

3. Phonon spectral function 

3.1 Normal state 

In our calculation, we typically set \g q \ 2 = 0.1, u> q = 
0.15 and consider the low temperature case T = 0.001, 
where A(T) ~ A(T = 0). We present the phonon spec- 



Dispersive Gap Mode of Phonons in Anisotropic Superconductors 



3 



tral function in the normal state A(T) = in Fig. 2(a). 
It shows S(q,co) along the ^-direction around Q. By 
the effect of Xnn(q,u) in the normal state, the phonon 
dispersion is shifted to lower energy than u) q , and it has 
a minimum at Q reflecting the nesting properties. The 
peak width of S(q,u>) comes from Imxnn(<7, The tem- 
perature dependence of the phonon dispersion is weak in 
this parameter case. The minimum energy of the disper- 
sion is cj = 0.088 at T = 0.05. It slightly decreases to 
lu = 0.086 at T = 0.001. 




Fig. 2. Phonon spectral function S(q = (q x ,0,0),ui) in the range 
0.4 < q x /n < 0.7 for the normal state A = (a) and the isotropic 
s-wave superconducting state with A = 0.025 (b), 0.035 (c) and 
0.050 (d). 



When ui q is further decreased, the phonon dispersion 
at q ~ Q approaches zero energy on lowering tempera- 
ture by the Kohn anomaly. This phonon softening case 
is considered in §4. 

3.2 Isotropic s-wave superconducting state 

We consider how S(q,u)) is changed in the supercon- 
ducting state. First, we investigate the isotropic s-wave 
case. We consider the A-dependence of the phonon prop- 
erties, where A = A(T = 0). Figures 2(b)-(d) show 
the change of S(q,u) with increasing A. When 2A ap- 
proaches the phonon dispersion from lower energy, the 
new peak of the gap mode appears below the phonon dis- 
persion, as shown in Fig. 2(b). The new peak is localized 
at q x <~ Q Xl when it appears. With further increasing 
2A, the gap mode extends to wider q x region, as shown 
in Fig. 2(c). The peak position of the gap mode shifts 
to higher energy when q x is away from Q. The original 
peak along the phonon dispersion is smeared above the 
gap mode. When 2A becomes larger compared with the 
phonon dispersion, the original peak at the phonon dis- 
persion vanishes, and the peak of the gap mode further 
grows up. To see this A-dependence clearly, we show 
the w-profile of S(q,oj) at q x = 0.55n in Fig. 3. The 
profiles at other q x also show the similar A-dependence, 
while characteristic A value is changed. It is also noted 
that the lower energy tail of the peak is suppressed by 
increasing A. 




CO 

Fig. 3. Profiles of the phonon spectral function S(q, u>) at q = 
tt(0.55, 0, 0) for A = 0, 0.02, 0.03, 0.04 and 0.05 in the isotropic 
s-wave case. 




A 0.05 



Fig. 4. A-dependence of the peak energy (upper panel) and the 
peak intensity (lower panel) for q x /n = 0.50 (o), 0.55 (•), 0.60 
(*). The points connected by solid lines correspond to the peak 
of the gap mode. Dashed lines are for the peak of the original 
phonon dispersion. We also show the line lo = 2A. 



In Fig. 4, we plot the A-dependence of the peak en- 
ergy and the peak intensity. We also show the cases 
q x = 0.507T and 0.607T in addition to the case q x = 0.557T 
in order to see the g x -depcndcnce. The peak intensity of 
the original phonon dispersion at each q x is monotoni- 
cally decreased with increasing A, while the peak energy 
at the phonon dispersion is not changed. When 2A ap- 
proaches the phonon dispersion, these original peaks are 
smeared, and the new peak of the gap mode appears at 
lower energy. It is seen as if the new peak first appears 
at finite lo. If we could exclude any contributions of peak 
broadening, we might observe new peak at lo ~ 2A until 
the lower A case near A <~ O. 11 ^ However, this new peak 
is not observed in the lower A case, since the peak in- 
tensity is rapidly decreased on lowering A and this peak 
is completely smeared out by the effects of peak broad- 
ening. Then, new peak appears at finite energy in our 
results. 

Compared with the cases q x — 0.57T and 0.6-7T, the new 
peak appears from smaller A at q x = 0.557T, since the 
phonon dispersion is located in lower energy at q x = 
0.557T. The increasing rate of the new peak's energy on 
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raising A is much weaker than the relation lo ~ 2A. At 
larger A, the peak energy is apparently smaller than 2A, 
and it seems to reduce to an ^-independent value which 
is smaller than the original phonon dispersion. With in- 
creasing A, the new peak becomes sharper and it has 
larger intensity. Along the dispersion for q x < Q x , both 
the original peak and the new peak are broad and over- 
lapped when new peak appears, as shown in Fig. 2(c). 
Because of the overlap, the energy of the new peak at 
q x = 0.57T in Fig. 4 is seen as if it appears at larger to 
and decreases with increasing A at first. 




q x /n 



Fig. 5. (/^-dependence of the crossover A, where peak intensity of 
the gap mode becomes larger than that of the original phonon 
peak, in the isotropic s-wave case (•), the s x 2_ y 2 case (o), and 
the s X y case (*). 

In Fig. 5, we plot the ^-dependence of the crossover 
A, where the peak intensity of the gap mode becomes 
larger than that of the original phonon peak. Above the 
line of the figure at each q x , the new peak of the gap 
mode is eminent. The crossover A is low near q x ~ Q x , 
and increases as q x is away from Q x . 

The A-dependence of the phonon profile in Fig. 3 
nicely reproduces the experimental data on YNi 2 B 2 C. 3 ~ 5 ) 
There, the A-dependence is the temperature dependence 
A(T) or the magnetic field dependence A(H). In the ex- 
perimental data, the energy of the new peak is close to 
2 A, but the T-dependence of the peak energy u> is weaker 
than the relation ui = 2A(T). This A-dependence of the 
peak energy is consistent with our results. The peak 
appears at the energy near 2A, when 2A approaches 
the phonon dispersion. When A further increases, the 
increase of the peak energy is weaker than the rela- 
tion to = 2A, as shown in Fig. 4. The experiment on 
LuNi 2 B 2 C reports that the phonon peak, which is shifted 
to lower energy, becomes sharp in the superconducting 
state. 7 ) It is consistent to our result. 

3.3 Anisotropic s-wave and al x 2_ y 2-wave superconduct- 
ing states 

In this subsection, we consider the anisotropic s-wave 
cases s x 2_ y 2 and s xy in order to clarify the effect of 
the superconducting gap suppression around the node 
points. While they have low energy excitations near 
nodes, we see that S(q, uj) shows qualitatively same be- 
havior as in the isotropic s-wave case presented in the 
previous subsection. That is, no drastic changes occur on 
phonon modes by their anisotropy. However, they show 
quantitatively different behaviors in the A-dependence 
and in the ^-dependence. We discuss these points in 



this subsection. 

The calculation of the d x 2_ y 2-wsve case gives quanti- 
tatively almost the same S(q,u) as that of the s x 2_ y 2 
case. This is because the nesting translation Q docs 
not change the sign of the pairing function <fid x 2_ y 2 (&) as 
shown in Fig. 1(c), and \4>d x2 _ y 2 (^)l nas the same gap 
amplitude as \<p s 2 _ 2 (fc)|. 




Fig. 6. Phonon spectral function S(q = (q x ,0,0),u>) in the range 
0.4 < q x /n < 0.7 for the anisotropic s-wave cases. In the s x 2_ y 2 
case, A = 0.020 (a), 0.030 (c) and 0.045 (e). In the s xy case, 
A = 0.025 (b), 0.035 (d) and 0.050 (f). The d x 2_ y 2-wmve case 
gives the same structure as that of the s x 2_ y 2 case. 

In Fig. 6, we present some typical behaviors of S(q, u) 
in the s x 2_ y 2 case and the s xy case. They show similar 
behavior as in the isotropic s-wave case in Fig. 2. When 
2A approaches the minimum of the phonon dispersion, 
a new peak of the gap mode appears at q = Q [(a) and 
(b)]. With increasing A, the new peak extends to the 
wide g^-range [(c) and (d)]. When A further increases, 
the original peak at the phonon dispersion is completely 
smeared out, and the peak of the gap mode becomes 
sharper [(e) and (f)]. The A-dependence of the peak 
energy and intensity is plotted in Fig. 7. The effect of 
the gap anisotropy appears at the onset A of the new 
peak and at the ^-dependence of the peak intensity. As 
for the onset A, we compare Figs. 4 and 7. In the 
case of the anisotropic pairings, the superconducting gap 
is enhanced or suppressed by the factor <f>(k) in A& = 
A(T)cj)(k). In the s x 2_ y 2 case, 0s x2 _ 2 (&) i s larger than 
1 along the flat Fermi surface < 8 < 30° as seen in Fig. 
1(b). At 9 = 0, 0s x 2_ !) 2 (k) has maximum 1.65. Then, we 
also plot the line u) = 2max(A^) = 3.3A in addition to 
the line uj = 2 A in Fig. 7. In Fig. 7(a), the onset A of 
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the new peak is lower than that of the isotropic s-wave 
case, but it is larger than the relation ui = 3.3A. That 
is, since the effect of (f) s 2 _ 2 (k) is averaged on the flat 
Fermi surface, the enhancement factor is larger than 1, 
but smaller than 1.65. Compared with the s x 2_ y 2 case, 
the onset A is larger in the s xy case, since the average of 
the factor Sxy (k) on the flat Fermi surface is smaller. 



0.1 



occurs. 20 -* The p x -wave case also shows the same behav- 
ior due to the sign change of </> Px (k) through the nesting 
translation. 
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Fig. 7. A-dcpcndcncc of the peak energy (upper panels) and the 
peak intensity (lower panels) for q x /ir = 0.50 (o), 0.55 (•), 0.60 
(*) in the s x 2_ y 2 case (a) and in the s xy -case (b). The points 
connected by solid lines correspond to the peak of the gap mode. 
Dashed lines are for the peak of the original phonon dispersion. 
We also show the lines ui = 2A and ui = 3.3A. 

As for the ^-dependence of the peak intensity, the 
intensity of the gap mode in Fig. 6 is enhanced for q x < 
Qx (qx > Qx) in the s x 2_ y 2 case (s xy case) compared 
with the isotropic s-wave case. The gap mode eminently 
extends toward q x < Q x (q x > Q x ) in Fig. 6(c) (Fig. 
6(d)). These characteristics are also recognized in Fig. 
7, where the peak intensity of the gap mode at q x = O.bir 
(0.67r) is largely enhanced, and in Fig. 5, where crossover 
A is lower for q x < Q x (q x > Q x ). This difference comes 
from the fc-dependence of <f>(k) along the perpendicular 
direction to the flat Fermi surface. We explain it by the 
fcc-dependence near the flat Fermi surface at k x = 0.2757T 
in Fig. 1(a). In the s x 2 _ y 2 case, the factor (j> Sx2 _ 2 C 2 ) is 
enhanced when k x is decreased because of the term cos k x 
in 0s x 2_ 2 Then, for the smaller q x (~ 2k x ) than the 
nesting vector Q x , the superconducting gap is enhanced, 
and the intensity of the gap mode becomes eminent. On 
the other hand, 4> Slcy (k) is enhanced when k x is increased 
because of the factor sinfc x in (j>s xy {k). Then, for larger 
q x (~ 2k x ) than Q x , the gap mode is eminent. 

3.4 d xy -wave and p x -wave superconducting states 

In this subsection, we examine the effect coming from 
the sign change of the pairing function by the nesting 
translation. The d xy -wa,ve case has the same supercon- 
ducting gap as that of the s xy case. But, the sign of the 
pairing function <pd xy (&) is changed by the nesting trans- 
lation Q as shown in Fig. 1(d). Then, the phonon be- 
havior in the d xy -~w&ve case is qualitatively different from 
the s X y case. It is because R,e{Afc-|_qA^} of \nn 

(q,uj) in 

cq. (2) becomes negative at q ~ Q when the sign change 




Fig. 8. Phonon spectral function S(q = (q x ,0,0),u>) in the range 
0.4 < q x /ir < 0.7 for the d^^-wave case (a) and the p x -wave case 
(b). A = 0.030. 

Figure 8 shows the typical behavior of S(q, ui) in the 
d^-wave (a) and the p^-wave (b) cases. There, the new 
peak of the gap mode does not appear. By the effect 
of A, the peak intensity of the original phonon mode is 
modified. But, the peak energy along the phonon dis- 
persion is almost unchanged. To see this A-dependence, 
we show the w-profile of S(q,u>) at q x = 0.557T for 
the dxy~ wave case in Fig. 9. With increasing A, the 
low energy tail is gradually suppressed. After 2A ex- 
ceeds the phonon dispersion, the peak intensity becomes 
sharper and the intensity is enhanced with increasing 
A. To see the ^-dependence of this behavior, we plot 
the A-dependence of the peak intensity at q x /ir = 0.50, 
0.55, 0.60 in Fig. 10. At low A, the peak intensity at 
q x = 0.507T (O.6O71-) is suppressed (enhanced) with in- 
creasing A, as also shown in Fig. 8. This g x -dependcncc 
comes from the fc^-dependence of the pairing functions 
\4>d xy {k)\ and \(j) Pic (k)\ near the flat Fermi surface at 
k x = 0.2757T. There, \(f>d x (k)\ and |0 Px (fe)| is enhanced 
with increasing k x because of the factor sin k x in the pair- 
ing functions. Then, for larger q x {^ 2k x ) than Q x , the 
superconducting gap is stronger, and the peak intensity 
is enhanced. 



500 



S 

5? 




Fig. 9. Profiles of the phonon spectral function S(q, u>) at q = 
tt(0.55, 0, 0) for A = 0, 0.03, 0.04 and 0.05 in the d^-wavc case. 



4. Suppression of phonon softening 

When ui q is enough low, the phonon dispersion ap- 
proaches zero energy at q = Q on lowering temperature 
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A 0.05 A 0.05 



(a) (b) 

Fig. 10. A-dependence of the peak intensity for q x /n = 0.50 (o), 
0.55 (•), 0.60 (*) in the d xy -wave case (a) and the p x -wave case 
(b). 



due to the Kohn anomaly. This phonon softening causes 
the second order transition to the CDW states. It is 
suggested that this softening is suppressed by the onset 
of superconductivity. 11 ^ To study this phenomenon, wc 
consider the phonon properties in the case ui q = 0.117. 
The other parameters are unchanged. 




Fig. 11. Phonon spectral function S(q = (q x , 0, 0), to) in the range 
0.4 < q x /w < 0.7 for w q = 0.117. In the normal state(A = 0), 
T = 0.050 (a) and 0.005 (b). In the superconducting state of the 
s^y-case (c) and the d xy -wave case (d), A = 0.05 and T = 0.005. 



Figure 11 shows S{q,u) in this case. At high temper- 
ature T = 0.050 (a), the phonon dispersion is located at 
finite energy. With decreasing temperature, the phonon 
dispersion is shifted to lower energy within the narrow 
region at q = Q. At low temperature T = 0.005 (b) in 
the normal state, the minimum of the dispersion touches 
zero energy, suggesting CDW transition. However, in the 
superconducting state, the phonon dispersion remains at 
finite energy, meaning the suppression of the CDW tran- 
sition. Figure 11(c) is for the s xy case. We obtain this 
type phonon behavior in the isotropic and anisotropic 
s-wave cases, and the rf a; 2_ y 2-wave case, i.e., in the case 
when the sign of the pairing function is not changed by 
the nesting translation. When the sign is changed by 
the nesting, such as the case d xy -wave and p^-wave cases, 
the phonon dispersion is further shifted to higher energy. 
We show S(q, uj) of the d xy -wa,ve case in Fig. 11(d). It is 
noted that these change of the phonon dispersion occurs 



only near q = Q. The dispersion is not largely changed 
far away from Q, e.g., at q x = OAtt or 0.7ir in Fig. 11. 

We plot the temperature dependence of the minimum 
energy of the phonon dispersion in Fig. 12, where the 
phonon dispersion is identified by the peak energy. In 
the normal state, the minimum energy ui monotonically 
decreases on lowering temperature. At T ~ 0.005, the 
CDW transition occurs by the softening u — > 0. On 
the other hand, the softening stops in the supercon- 
ducting state. Then, the second order transition to the 
CDW state does not occur. The minimum w shows weak 
temperature dependence below T c in the s-wave cases 
and the d x 2_ y 2-v?8cve case. However, the minimum uj is 
largely enhanced just below T c and it saturates at low 
temperature in the d xy -wa.ve and the p^-wave cases, i.e., 
in the case when the pairing function changes the sign 
through the nesting translation. In Fig. 12, T c is differ- 
ent due to the factor cj>(k), while we set A(T = 0) = 0.05 
in each pairing case. 




Fig. 12. Temperature dependence of the minimum energy of the 
phonon dispersion for ui q = 0.117 in the normal state (A = 
0) and the superconducting state with A(T = 0) = 0.05. In 
the superconducting state, we present the isotropic s-wave, the 
anisotropic s-wave s x 2_ y 2 and s xy , the d xy -wecve and the p x - 
wave cases. The points for the d x 2_ y 2-wave is the same as the 
s x 2_ y 2 case. 

The amplitude of the superconducting gap is also con- 
trolled by the external magnetic field. In the case when 
the CDW transition is suppressed by superconductivity, 
the CDW state is recovered if we destroy the supercon- 
ductivity by applying the magnetic field. This predic- 
tion is testable for appropriate superconductors, possibly 
YNi 2 B 2 C. 

5. Summary and discussions 

By extending the previous theories 8-11 ) on the appear- 
ance of the new phonon mode below T c for isotropic s- 
wave superconductors, we have investigated the detailed 
features of the dispersion relation for this new phonon 
gap mode in light of various possible anisotropic pair- 
ings. It is demonstrated that observation of this dis- 
persive gap mode can be a useful spectroscopic tool for 
sorting out some of the possible pairing states. With 
increasing the amplitude of the superconducting gap A, 
this new peak of the gap mode becomes sharp and the 
original phonon peak is smeared. We have estimated 
the A-dcpcndcncc of the dispersion and the intensity of 
the peak. The effect of the gap anisotropy is classified 
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to two cases; (i) the case when the sign of the pairing 
function is not changed through the nesting translation, 
and (ii) the case when the sign is changed. In the former 
case (i), the gap anisotropy modifies the phonon spectral 
function quantitatively. The anisotropy effects appear at 
the onset A of the gap mode and at the ^-dependence 
of the peak intensity reflecting local amplitude of the 
superconducting gap at the Fermi surface region related 
to the nesting translation. In the latter case (ii), the 
new peak of the gap mode does not appear. There, for 
larger A, the original phonon peak becomes sharp with- 
out large shift of the peak energy. While we obtain these 
properties by analyzing example models for the Fermi 
surface and pairing functions, we can expect that these 
qualitative results coming from the superconducting gap 
anisotropy are independent of the model, because of the 
generality of their origin. 

In the experimental data on YNi 2 B 2 C, the eminent 
dispersive new peak of the gap mode is observed around 
q = Q. 3 ~ 5 ) This behavior of the gap mode is qualita- 
tively consistent to our numerical results. The phonon 
peak on LuNi 2 B 2 C is shifted to lower energy in the su- 
perconducting state. 7 -' Then, in these materials, we can 
exclude the pairing symmetry of the case (ii), for exam- 
ple, (f^-wave. That is, we can conclude that the pairing 
function has the same sign at the flat Fermi surface re- 
lated to the nesting. 

Further identification of the pairing symmetry among 
isotropic s-wave, anisotropic s-wave and d x 2_ y 2-wa,ve is 
difficult, because it needs quantitative considerations. 
But, we give some discussions here. One of the key points 
to identify the pairing symmetry is the ^-dependence of 
the peak intensity. In the data on YNi 2 B 2 C, the peak 
intensity is larger for q x < Q x . 3 ^ Following our anal- 
ysis, this suggests that the gap anisotropy of the type 
cos fez — cos k y \ is preferable. The second point is that 
the superconducting gap largely opens at the flat Fermi 
surface related to the nesting, when the eminent gap 
mode appears. The contribution of the superconduct- 
ing gap affects on the phonon as an average along the 
flat Fermi surface. In our calculation, the average of 
(p Sx (fe) is smaller than that of (f> s 2 _ 2 (k). The former 
has a node at the flat part of the Fermi surface. Since 
the flat part is long in our model of the Fermi surface, 
the both averages do not show large difference. How- 



ever, if the flat part is shorter, the average of the gap is 
smaller in the case when the node of the gap anisotropy 
is located at the fiat part. In this case, we cannot expect 
the eminent gap mode. To obtain decisive conclusion on 
these above-mentioned points, we need to establish the 
knowledge about the detailed shape of the Fermi surface 
and the nesting properties. 

Note added: Recently, a similar neutron experiment 
has been reported on ErNi 2 B 2 C by Kawano-Furukawa 
et al, showing the gap mode below T c . 21 ) 
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